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^> ' Abstract 

m 

. Despite a fainter Sun, the surface of the early Earth was mostly ice-free. Pro- 

I posed solutions to this so-called "faint young Sun problem" have usually involved 

^0 . higher amounts of greenhouse gases than present in the modern-day atmosphere. 

^ I However, geological evidence seemed to indicate that the atmospheric CO2 concen- 

^ ■ trations during the Archaean and Proterozoic were far too low to keep the surface 

I from freezing. With a radiative-convective model including new, updated thermal 

absorption coefficients, we found that the amount of CO2 necessary to obtain 273 K 
at the surface is reduced up to an order of magnitude compared to previous studies. 
For the late Archaean and early Proterozoic period of the Earth, we calculate that 
CO2 partial pressures of only about 2.9 mb are required to keep its surface from 
freezing which is compatible with the amount inferred from sediment studies. This 
conclusion was not significantly changed when we varied model parameters such 
as relative humidity or surface albedo, obtaining CO2 partial pressures for the late 
Archaean between 1.5 and 5.5 mb. Thus, the contradiction between sediment data 
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and model results disappears for the late Archaean and early Proterozoic. 
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1 1 Introduction 



2 Geological evidence has shown tha t liquid water was present on the Eart 



face earlier than 3.7 Gy ago (e.g., 



h s sur- 



Mojzsis et al. 



1996 



Rosing and Frei 



20M) 



which implies average temperatures on the surface above 273 K. Some authors 



have even argued 



or a h ot Archaean climate (T >340 K), based on oxygen 



( iKnauth and Lowd . l2003l ) and silicon ( IRobert and Chauss 



analys is of seawater cherts. How e ver, a s pointed out by, e.g. 



donl. 



2nn6h isotope 



Kasting and Howard 



(120061 ) and [Shields and KastingI (|2007| ). these isotopic signature s chang es might 



not only be caused by temperature effects. iSleep and Hesslerl (120061 ). for ex- 

10 ample, deduced a more moderate surface temperature below 300 K, based on 

11 quartz weathering records in paleosols. Nevertheless, it is generally accepted 

12 that the Earth has been ice-free throughout most of its history. 

13 Observations of several solar-type stars of different ages and virtually all stan- 

14 dard models of the solar interior have shown that the total solar lumino sity has 



increased since the ZAMS 



1981 



Zero Age Main Sequence) by about 30% (iGough 



Caldeira and Kasting] . Il992l ). Had the composition of the Earth's atmo- 



17 sphere been the same then as today, the reduced sol ar flux would have resulte d 



1972) 



18 in surface temperatures below 273 K prior to 2.0 Gy (ISagan and MuUenl . 

19 This apparent contradiction between solar evolution models, climatic simula- 

20 tions and geological evidence for liquid water and moderately warm tempera- 

21 tures on Earth has been termed the "faint young Sun problem". 



22 Numerous studies have attempted to solve this problem. For example. 



Minton and Malhotra 
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23 (120071) explor ed the hot early Sun scenario for a non-standard solar evolution. 



Shavivl (120031 ) showed that moderate greenhouse warming in combination with 



25 the i nfluence of sola r wind and cosmic rays on climate could resolve the prob- 



lem. 



Jenkins! (120001 ) assumed high obliquities in a General Circulation Model 



27 to account for high Archaean temperatures. 



However, the most accepted scenario involves a much enhanced greenhouse 
effect (GHE) on the early Earth compared to modern Earth. Today, the GHE 
produces around 30 K of warming, raising the mean surface temperature of the 
Earth to abou t 288 K. Incr eased abu ndances of greenhou se gases such as car- 



32 bon dioxide (iKastingl. 



2008! ) or ammonia (jSagan and Mullen 



19871). methane ( Pavlov et al 



1972 



200o[). ethane (iHaqq-Misra et al. 



Sagan and Chybal . 119971 ) will strengthen 



the GHE, hence potentially resolve the problem. However, all of these studies 
faced some form of contradictions or large uncertainties, either from geological 
data on atmospheric conditions or from atmospheric modeling. The formation 
and des truction of ammonia is highly dependent on UV levels in the atmo- 



sphere ( ISagan and Chyba 



1997 



Pavlov et al. 



200ll ). The hydrocarbon haze 



39 necessary to allow higher hydrocarbons to accumulate in t 



pends critically on the CO2/CH4 ratio (jPavlov et al. 



le atmosphere de- 



20031 ). The high values 



of methane required to heat the surface of the early Earth depend on estimates 
of the early biosphere and volcanic activity, which is n ot well d e termi ned. Past 
CO2 concentrations required by atmospheric models (IKasting) . 119871 ) to reach 
surface temperatures ab ove 273 K are in conflict with inferre d concentrations 



from the sediment data (IHessler et al. 



2004 



Rye et al. 



19951). 



46 In this work, the role of CO2 in warming the early Earth is reconsidered. We 

47 used a one-dimensional radiative-convective model, including updated absorp- 

48 tion coefficients in its radiation scheme for thermal radiation. The model was 

49 applied to the atmosphere of the early Earth in order to investigate the effect 

50 of enhanced carbon dioxide on its climate. Additionally, we investigated the 

51 effect of two important parameters, namely the surface albedo and the relative 
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52 humidity, upon the resulting surface temperature. 

53 Our results imply that the amount of CO2 needed to warm the surface of early 

54 Earth might have been over-estimated by previous studies. Furthermore, the 

55 results show that the contradiction between modelled CO2 concentrations and 

56 measured values might disappear by the end of the Archaean. 

57 Section [2] describes the model and section [3] the runs to validate the new radia- 

58 tion scheme. The runs performed for this work are explained and summarized 

59 in section HI In section the results are presented and discussed. Section [H] 

60 gives the summary of the results. 



61 2 Atmospheric Model 



62 We used a one-dimensional 

63 part of the model used by 



radiative-convective model based on the climate 



Segura et al.l (120031 ) and 



Grenfell et al. 



fl2007a 



64 Our model differs in upgrades of the radiation scheme to calculate the thermal 

65 emission in the atmosphere. The model calculates globally, diurnally-averaged 

66 atmospheric temperature and water profiles for cloud-free conditions. We will 

67 first state some basic characteristics of the model (12.11) . Then we will describe 

68 the calculation of the energy transport via radiative transfer (solar and thermal 

69 fluxes) and convection (12.21) to obtain the atmospheric temperature profile. 

70 Thereafter, a description of the determination of the water profile is given 

71 (12. Sp . Finally, the model input parameters are summarized (12.41) . 



72 2.1 Basic model description 



73 The model determines the temperature profile by assuming two dominant 

74 mechanisms of energy transport, i.e. radiative transfer and convection. The 

75 convective lapse rate is assumed to be adiabatic. The radiative lapse rate is 
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76 calculated from contributions of both solar and thermal radiation, including 

77 Rayleigh scattering for solar radiation and continuum absorption in the ther- 

78 mal region. Table [T] summarizes the contributions of the different atmospheric 

79 species to the calculation of the temperature profile. 

80 The species considered in the model are molecular nitrogen, water, molecular 

81 oxygen, argon, carbon dioxide and carbon monoxide. For example, a typi- 

82 cal early Earth run considered 0.77 bar of nitrogen and several (2.9 - 57.2) 

83 mb of carbon dioxide in addition to water in varying concentrations (0.5 — 1 

84 % at the surface). Molecular nitrogen is an effective Rayleigh scatterer (al- 

85 though not as effective as carbon dioxide) and as a main constituent of the 

86 atmosphere also contributes to the heat capacity. Water is not considered to 

87 be an important Rayleigh scatterer, but it is relevant for the other radiative 

88 processes. Also, water influences the adiabatic lapse rate because it readily 

89 condenses in the troposphere. However, due to small mixing ratios, especially 

90 in the stratosphere, water vapour does not contribute to the heat capacity of 

91 the atmosphere. Molecular oxygen and argon do contribute to the heat capac- 

92 ity of the atmosphere, and molecular oxygen additionally contributes to the 

93 Rayleigh scattering coefficient. Carbon dioxide contributes to all relevant ra- 

94 diative mechanisms (molecular absorption of solar and thermal radiation, con- 

95 tinuum absorption, Rayleigh scattering), but not to the adiabatic lapse rate 

96 because it does not condense under conditions described in this paper. Carbon 

97 monoxide is an important absorber species in some mid-infrared windows and 

98 contributes to the heat capacity. 

99 The model assumes the hydrostatic relation between pressure p and density 

100 g throughout the plane-parallel atmosphere. On the 52 model layers, a loga- 

101 rithmic pressure grid is calculated. Specified pressure levels at the planetary 

102 surface (e.g., 1 bar for the standard Earth case) and the upper model lid (at 

103 6.6-10"^ bar) determine the altitude range which, for modern Earth conditions, 

104 extends to 65-70 km, i.e. the lower- to mid mesosphere. For all gases except 
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105 water, the ideal gas law is taken as the equation of state (see [2231 for water). 



106 The effect of c 



Pavlov et al. 



approach of 



ouds is difficult to incorpo rate into ID models (see, for example. 



2000 



Segura 



Segura et al. 



eta 



20031 ). In the present model, following the 



(120031 ). clouds are implicitly included by adjusting 
109 the surface albedo Agurf such that under modern Earth control conditions the 
no model calculates a mean surface temperature of 288 K. The required value for 

111 v4surf is about 0.21, whereas the actual global value for Earth is approximately 

112 0.15. This can be interpreted as a ground cloud layer instead of tropospheric or 



113 stratospheric clouds. The m odel uses a time-ste 



114 to the steady-state solution ljPavlov et al.l . 120001 ). 



Dping algorithm to convergence 



115 2.2 Temperature profile 



116 During each time step, the temperature profile is calculated from the radia- 

117 tive equilibrium condition. The temperature T on an atmosph eric level z is 



118 determined by the following equation of energy conservation (jPavlov et al 



2000): 



dt ^ ' Cp{T, z) dp{z) ^ ' 

120 where dt is the time step in the model, Cp the heat capacity, F the total net 

121 radiative flux and p the pressure of the level. The radiative flux F is the sum of 

122 thermal planetary and atmospheric emission, Fthermai, and the solar radiative 

123 input, Fsoiar, into the atmosphere: 



F{z) = FthcrmaK^;) + -^solar ( ^ ) (2) 

124 These fluxes are calculated separately by two numerical schemes which solve 

125 the monochromatic radiative transfer equation (RTE) for the spectral intensity 

126 lu in the respective spectral domain (i.e., near-UV to near-IR for solar flux. 
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127 near-IR to far-IR for thermal flux): 



^-T- = h- (3) 

128 where is the source function (either the incident solar radiation or the ther- 

129 mal blackbody emission of the atmospheric layers and the planetary surface), 

130 dr^ the optical depth and /i = cos(^) the cosine of the polar angle. The optical 

131 depth is defined as usual by 



dr = -(/Cy + s,,)dz (4) 

132 where and represent the absorption coefficient and the scattering coef- 

133 ficient respectively. The absorption coefficient for a gas mixture is calculated 

134 from the individual absorption coefficients of the gas species i: 



l^u = Y1 ^'^^i = H ^abs,i ■ Ni (5) 

i i 

135 where is the number density and (Jabs,i the molecular absorption cross section 

136 of the gas species i. When no scattering occurs (i.e., Sy = 0), eq. H] can be 

137 written in terms of the column density Wi of the gas species i as: 

r = X! ^abs,i • Wi (6) 

i 

138 The absorption cross section is defined by: 



(Tabs (i^, p,T) = Y, Sj (T) ■ gj{u, T, p) (7) 
j 

139 Here, Sj{T) is the temperature-dependent line strength of a particular spec- 

140 tral line j and gj{h',T,p) the temperature- and pressure-dependent line shape 
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141 function of the same line. 



142 For the scattering coefficient, an analogous equation is valid: 



i i 

143 Here, CTyA^) is a scattering cross section of type y. In the present model, 

144 Rayleigh scattering is considered. 

145 The solution of eq. (j4]), i.e. the calculation of optical depths, is one of the key 

146 elements in radiative-convective models. This solution, of course, depends on 

147 the accurate calculation of the absorption cross sections. 

148 From eq. ([3]), the necessary fluxes for eq. ([2]) (i.e., the thermal and the solar 

149 flux) are obtained by an angular integration of the (monochromatic) intensity: 



dujjily (9) 
150 and a frequency integration of the (monochromatic) flux: 



rv2 

F,,^,, = / F,du (10) 

151 Each one of these integrations is performed independently for the two compo- 

152 nents of the total flux. 



153 2.2.1 Solar radiation 



154 The solar radiatio n module whi c h cal c ulates Ff^rA^Az ) for e g. ([2]) has already 



155 been u sed by, e.g. 



Pav 



156 (l2007al ) and is based on 



ov et al 



fl2nnnh: 



Kasting et al 



Seeura et 



fll984b[ ) and 



al 



fl2nn3h or 



Grenfell et al 



KastineJ fllQSSh . The module 



157 considers a spectral range from 0.26 to 4.5 /im in 38 intervals. It evaluates the 
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158 solar incident radiation at a fixed daytime average solar zenith angle of 60°. 

159 Contributions to the optical depth come from gaseous absorption by water 

160 and carbon dioxide (i.e., in eq. (jlj)) and from Rayleigh scattering by carbon 

161 dioxide, molecular nitrogen and molecular oxygen (i.e., in eq. (jlj)). 



162 Absorption cross sectio ns CTahs for the solar c ode were obtained from the HI- 



163 TRAN 1992 database (jRothman et al. 



1992j). Rayleigh scat t ering cross sec- 



Vardavas and Carver 



(jl9M). The fre- 



164 tions a ray are parameterized following 

165 quency integration (see also eq. ( fTOl) ) of the RTE for Fgoiar in each of the 38 

166 spect ral intervals is parameterized by a four-term correlated-k exponential sum 



167 (e.g. 



Wiscombe and Evans 



19771 ). The following angular integration (eq. 



168 is performed by u sing a quadrature 5-2-stream approximation code based on 



Toon et al. 



(jl989l ). The resulting fluxes from each spectral interval are added 

170 up to yield the total solar flux Fsoiar(-2) at an atmospheric level z. This flux is 

171 further multiplied by a factor of 0.5 to account for diurnal variation. 



172 2.2.2 Thermal molecular absorption 



The thermal (planetary) radiation module for Fthcrmai(-2) in eq. (j2J) consid- 
ers a spectral range from 1 to 500 fim in 25 intervals. Our new thermal 
module is called MRAC (Modified RRTM for Application in C02-dominated 



176 Atmospheres) and is based on the radiation s cheme RRTIV 



tive Transfer Model). RRTM was developed by 



Mlawer et al. 



been used by numerou s other modeling studies (e.g. 



Rapi d Radia- 



Grenfell et al 



2007a 



Segura et al. 



f|l997h and has 



2003 



2005 



b|). The need for a new radiation model comes from the 



fact that RRTM was specifically designed for conditions of modern Earth, i.e. 
it is not adaptable to studies of atmospheres which greatly differ from mod- 
ern atmospheric conditions (in terms of atmospheric composition, temperature 
structure, pressure, etc.). MRAC is easily adaptable to varying conditions, as 
is described below. 



185 MRAC uses the correlated-A; approach (e.g., 



Goody et al. 



1989 



Lacis and Oinas 
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186 119911 : iColaprete and Toonl . 120031 ) for the frequency integration of the RTE in 

1B7 the thermal range, as does RRTM. The planetary surface (bottom layer of the 

1B8 model atmosphere) and the atmospheric layers are taken as blackbody emit- 

1B9 ters, according to their respective temperatures. The thermal surface emissivity 

190 is set to unity. The absorber species considered in MRAC are water, carbon 

191 dioxide and carbon monoxide. The angul ar integration (eg. (p) ) is performed 



192 using the diffusivity approximation (as in 



Mlawer et al. 



1993). 



193 The correlated-k method 



194 Basically, this method transfers the frequency integration (FI) of the RTE from 

195 frequency space to a probability space g. For the absorption cross section 
CTabs of eq. ([7]) in an interval [ui, z/2], a probabilistic density distribution f (crabs) 

197 (i.e., probability of occurence for a particular value of o"abs) with the following 

198 normalization condition can be defined: 



poo 

/ /(0-abs)<i0-abs = 1 (H) 

Jo 

199 To follow conventional nomenclature in the literature regarding /c-distributions, 
CTabs is hereafter referred to as k. The function /(cXabs) = f{k) is called the k- 

201 distribution. From the /c-distribution, a cumulative /c-distribution g{k) can be 

202 defined by 



g{k) = I f{k')dk' (12) 
Jo 

203 The cumulative ^-distribution is a strictly monotonic function and may thus 

204 be inverted from g{k) to yield k{g). This mapping of the frequency informa- 

205 tion {k{v)) into a single probability variable {k{g)) can be done because it is 

206 irrelevant at which position of the spectral interval a particular value of the 

207 absorption cross section k occurs. Performing a variable substitution in eq. 

208 ( fTOl) then leads to the following equation: 
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2^2 



F{v)du = / F{g)dg 

VI Jo 



(13) 



209 The goal of the cumulative-fc approach is to reduce the number of radiative 

210 transfer operations drastically while keeping the accuracy of line-by-line mod- 

211 els. This can be a c hieved with very few numbers of points in g space (see 



Goody et al. 



1989 



West et al. 



1990). The g integration in eq. (1131) is per- 



213 formed in MRAC by G aussian quadrature u sing 16 intervals in g space (the 



214 same as used in RRTM, 



Mlawer et al. 



19971 ). 



215 The extension of this exact method from homogen eous to inhornogene ous at 



Mlawer et al. 



19971). Each 



216 mospheres is called the correlated-A; method (e.g., 

217 g interval is treated as if it were a monochromatic frequency interval, i.e. the 

218 method uses the same subset of g space for all layers throughout the atmo- 

219 sphere. This implies a full frequency correlation of a specific subset of g space 

220 for all atmospheric layer s . The re are conditions under which this approach is 



221 exact (iGoody and Yung 



1989I ). but usually these do not hold. However, the 



222 nume rical error of the correlated-A; method is generally small (IMlawer et al. 



19971). 



Creating the new radiation scheme 



225 MRAC was originally designed to simulate atmospheres of a wide range of pos- 

226 sible terrestrial planets other than modern Earth, as stated above. Therefore, a 

227 new temperature-pressure (T-p) grid to incorporate the aforementioned {T,p)- 

228 dependence of the absorption cross sections (see eq. [7]) has been introduced. 

229 In RRTM, k values are tabulated for every spectral band and every point in g 

230 space for 59 pressure levels and the associated Mid-Latitude-Summer (MLS) 

231 standard Earth temperature valu es as well as tempera ture values Tmls± 15 K 



232 and Tmls± 30 K, as described in 



Mlawer et al. 



(119971 ). The T-p grid of RRTM 



233 is thus more or less fixed to modern Earth conditions. For MRAC, we used 

234 8 pressure levels, ranging equidistantly in log p from 10~^ to 100 bar and 9 



11 



235 temperature points, 6 in 50 K steps from 150 K to 400 K and three additional 

236 points for 500, 600, 700 K respectively. Tests showed that this grid allows for 

237 an interpolation accuracy of usually better than 2-3%. Furtherr nore, the num- 



238 ber of T-p points is consistent w ith previous modeling studies (IKasting et al. 



1993 



Colaprete and Toonl . |2003| ). 



240 Figure [T] shows the range of the tabulated /c-values for both RRTM and MRAC, 

241 i.e. the interpolation regime of the two radiative schemes. It demonstrates that 

242 RRTM can only be applied to a much narrower range of atmospheric conditions 

243 in comparison to MRAC. 

244 The necessary re-calculation of the /c-distributions for each of the gases in- 

245 eluded in MRAC (water, carbon dioxide, carbon monoxide) proceeded in three 

246 steps: 

247 First, the absorption cross sections for every species were calculated for each 

248 T-p grid point in each of the spectral intervals where the respective species is ac- 

249 tive. The line shape cut-off was set to 10 cm~^ from the frequency i.e. the sum 

250 over j in eq. (JTj) contains contributions from all lines within ± 10 cm'^. For the 

251 line shape fi^j (z^, T,p), a. Voigt profile was assu med. The required line p arameters 

252 were taken from the HiTemp 1995 database (IRothman et al.l . Il995l ). The for- 

253 eign broadening parameters in HiTemp are given for air, i.e. an oxygen-nitrogen 

254 mixt ure as a background atmo s phere . However, as reported by several authors 



255 (e.g. 



Brown et al. 



2005 



Toth 



2000), the foreign broadening parameters vary 



by significant amounts when different broadening gases are considered. Thus, 
for each type of background atmosphere, a new set of ^-distributions must be 
generated (e.g., low-oxygen atmospheres, C02-dominated atmospheres, inter- 
mediate N2-02-atmospheres). These line parameters we re then used as input 
to a line-by-line radiative transfer model called Mirart (jSchreier and Bottgerl . 



20031 ). Mirart produced the actual absorption cross sections with a spectral 



262 resolution of 10^ equidistant points per spectral interval. 



263 Second, the /c-distributions f{k) were calculated from the absorption cross 
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264 sections. From the k distribution f{k), the cumulative /c-distribution g{k) was 

265 then obtained. 



266 Third, representative k values were calculated for each of the q subinte rvals. In 



Mlawer et al 



fll997h . i.e. for 



267 this step, our algorithm followed the approach of 

268 each of the 16 Gaussian subintervals in g space, an arithmetic mean absorption 

269 cross section was calculated. 



270 MRAC also implements a so-called binary species parameter r] for transmit- 



271 tance calculations. This is e mployed in intervals 



272 species (for more details, see 



Mlawer et al. 



la ving two important a 



1997 



3Sorber 



or 



Colaprete and Toon 



20031 ): 



273 Here, Ci^2 are the concentrations of the two gases (in the case of MRAC, wa- 

274 ter and carbon dioxide) and r is some specified reference ratio (mean modern 

275 Earth tropospheric values). Carbon monoxide, although present in the model 

276 in six spectral intervals, is not considered to be part of the binary species 

277 parameter. This is partly due to the expected low concentrations in the sim- 

278 ulated atmospheres (typical theoretical and measured values for early Earth 

279 and Mars are below 10~^ volume mixing ratio), partly due to the expected low 

280 temperatures (below 300 K), which means that the strong CO fundamentals 

281 around 4 /im are completely outside the relevant Planck emission windows. 

282 Consequently, carbon monoxide has only a reduced impact on the radiation 

283 budget, compared to water and carbon dioxide, ^-distributions are calculated 

284 in MRAC for 5 different values of rj, ranging equidistantly from (CO2 only) 

285 to 1 (H2O only). 

286 Another improvement in MRAC, compared to RRTM, is the treatment of the 

287 Planck function in each band. The fraction of thermal radiance associated with 



288 a subset in g space is calculated from eq. (11) in lMlawer et al.l (119971 ): 
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(15) 



289 Here, Bg is the average Planck function of the frequencies in the subset of g 

290 space, Wg the Gaussian weight of the g interval and Bg the average Planck 

291 function of the whole spectral band (i.e., the whole g space). As temperature, 

292 pressure and species concentrations vary, the different g subsets will correspond 

293 to different frequencies, and as such the value of Bg, thus fg, will vary. This 

294 has been taken into account while constructing the /c-distributions. 



Mlawer et al.l (Il997l ) tabulated values of fg for every value of the binary species 

296 parameter and for two atmospheric reference levels, one each in the troposphere 

297 and the stratosphere. In MRAC, values of fg were tabulated for three temper- 

298 atures and two pressure levels as well as for the values of the binary species 

299 parameter. The most important factor for fg is the binary species parameter r], 

300 whereas the variation with temperature and pressure is rather small, although 

301 not negligible. Therefore these 6 T-p points are regarded as to be sufficient. 



302 A further difference between RRTM and MRAC is the distinction between 

303 troposphere and stratosphere. In the troposphere, RRTM chan ges maior anc 



304 mino r absorbers in some of the spectral intervals (see table 1 in 



Mlawer et al. 



19971 ). In some spectral bands, no absorption is considered in the stratosphere. 



306 in others, the number of key species is reduced. This is done because on Earth, 

307 the chemical and physical regimes are quite different in the troposphere com- 

308 pared to those in the stratosphere. However, as this is mostly due to Earth- 

309 specific conditions (e.g., the cold trap, tropopause and temperature inversion 

310 all approximately occur at the same altitude), this distinction was not incor- 

311 porated into MRAC. 
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312 2.2.3 Thermal continuum absorption 



313 Based o n ap proximation formulat i ons us ed by iKasting et al.l (1l984bl ) , iKasting et al. 



1984al ) and IColaprete and Toon! (120031 ). additional CO2 and H2O continuum 



315 absorption in the thermal region is considere d. In contrast, the R RTM scheme 

316 only considers water continuum absorption (IMlawer et al.l . 119971 ). 



317 Equation ( U61) shows the approximation for the optical depth rcont,co2 ^'^^ 

318 to CO^ continuum ab sorption. The corresponding parameters are taken from 



Kasting et al. 



fll984bh . 



^cont.COa = ^iW ■ Pe 



T 



(16) 



320 In this equation, Ci a frequency-dependent adjustment to the path length, 

321 W the column amount of CO2, pe = + 0.3 ■ CCO2) ' P ip layer pressure, 

322 Cco2 concentration) an effective CO2 broadening pressure and To=300 K is a 

323 reference temperature. 



324 The optical depth rcont,H20 due to water continuu m absorption in the w indow 

325 region (8-12 //m) is calculated from the equation (IKasting et al.l . Il984al ) 



7 w 

7"cont,H20 — iT-ii ■ P ■ -TTT 

Wt 



(17) 



326 where hn = hfhy incorporates the frequency and temperature dependence, p is 

327 the pressure, Wt the total and Ww the water col umn of the laye r , hy is e valuated 



328 at the high frequency interval boundary, as in 



Kasting et al. 



329 the 

330 and 



bllowing approxima t ions 



Colaprete and ToonI (120031 ): 



br ht and h^, based on 



1984a[). We use 



Kasting et al 



(Il984a[ ) 
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K = 1.25 • 10-2^ + 1.67 • 10-1^ • (^g) 

331 Both the water and the carbon dioxide continuum absorption are considered to 

332 be approximately monochromatic over a specific spectral interval, hence their 

333 contribution to the overall absorption coefficient (see eq. (jl])) is added as a 

334 constant term. 



335 2.2.4 Convective adjustment 

336 Convective adjustment to the lapse rate is performed whenever the calculated 

337 radiative lapse rate Vrad^" exceeds the adiabatic value Vad^ (Schwarzschild 

338 criterion) : 



(20) 



339 The adiabatic lapse rate is calculated as a standard dry adiabat in the strato- 

340 sphere. In the troposphere, a wet H2O adiabatic lapse rate is assumed. Below 

341 273 K, the Clausius-Clapeyron-equation '^Ji^l^^^ = {R universal gas con- 

342 stant, m mass, L latent heat release per mass) for the saturatio n vapor p r essure 



343 curve is applied. Between 273 and 647 K, a formulation by 

344 is taken. 



Ingersolll ( 119691 ) 



345 2.3 Atmospheric water profile 

346 In every time step, the water vapor profile is re-calculated according to the 

347 new temperature profile. 

348 In the troposphere, water vapor concentrations Cn^o are calculated from a 

349 fixed relative humidity distribution RH: 
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p{z) 



RHiz) 



(21) 



350 where Ps&t is the saturation vapor pressure of water at the given temperature T 

351 and p the atmospheric pressu re at level z. The default r elativ e humidity profile 

352 RH follows the approach of iManabe and Wetheraldl (119671 ). with a relative 

353 humidity Rg of 80% at the surface. 



RHiz) = R, ■ ^2^^ 



0.02 



0.98 



(22) 



354 Above the cold trap, water vapor is treated as a non-condensable gas, and its 

355 concentration is fixed at the cold trap value. 



356 2.4 Boundary conditions, initial values and parameters 



357 Since eq. ([T]) is a first order differential equation for the temperature, a start- 

358 ing temperature profile must be provided. In addition, a boundary condition 

359 for the radiative flux must be specified. To obtain unique equilibrium solu- 

360 tions, parameters must also be provided for the model. These include pressure 

361 parameters for the planetary top-of-atmosphere (TO A) pressure po, gas con- 

362 centrations, surface albedo or solar zenith angle, for example. 

363 Table [2] summarizes the boundary conditions, initial values and parameters. 



364 3 Validation of the new radiation scheme 



365 3.1 General remarks 



366 MRAC has been tested in two different ways. 
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367 • Case 1: /c-distributions 

368 The calculated /c-distributions, i.e. the model input data, have been val- 

369 idated against published values to show that the algorithm creating the 

370 /^-distributions works correctly. 

371 • Case 2: Earth atmosphere temperature profiles 

372 Temperature profiles of an Earth-like test atmospheres (composition: N2 

373 0.77, O2 0.21, Ar 0.01, CO2 3.55 ■ 10"^) calculated with MRAC and RRTM 

374 have been compared. This was done since RRTM has been extensively 

375 yalidated both agai nst line-by-line codes and atmospheric measurements 



376 ( iMlawer et al. 



19971 ) under modern Earth conditions. 
Our test atmosphere has a composition close to the present day atmo- 
sphere. However, it lacks radiative trace gases such as nitrous oxide, ozone 
and methane, as these gases cannot be handled by MRAC yet (see above). 
Note that due to the lack of ozone in our test atmosphere, we do not expect 
a large stratospheric temperature maximum as is observed in the present 
Earth atmosphere because this maximum is almost entirely due to the ab- 
sorption of solar radiation by ozone. 

We additionally performed test runs on a second test atmosphere (not 
shown) which differs from the first one by its CO2 content (100-fold increase). 
This 10 0-fold increase in C O2 represents the current limit for the RRTM 



scheme ( ISegura et al. 



2003 



, Eli Mlawer, priv. comm. 



388 These two validation approaches are discussed below. 



389 3.2 k- distributions 



390 Figures [2] and [3] compare our calculated /c-distributions (dotted lines) with 

391 previously published values (plain lines) fo r different water a nd ca rbon dioxide 



392 bands. Published values were taken from 

393 [2]) and 



Mlawer et al. 



Lacis and OinasI fll99lh (H2O, Fig. 



(119971 ) (CO2, Fig. ED. Figures [2] and [3] indicate quite good 



394 agreement with the published values. 
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395 3.3 Earth temperature profiles 

396 The calculated temperature profiles for the test atmosphere are shown in Fig. 

397 H 

398 Figure H] implies some differences in the middle to upper stratosphere (2-6 K) 

399 and small deviations (<C IK) below about 20 km. For the test atmosphere 

400 with a 100-fold increase in CO2 (not shown), the stratospheric differences are 

401 even larger (up to 10 K). We interpret these differences in the temperature 

402 profiles as follows: 

403 Firstly, as stated above in section [2.2.21 MRAC does not differentiate between 

404 troposphere and stratosphere, as is the case for RRTM. That means, spectral 

405 bands where H2O or CO2 absorb only weak are not considered for optical depth 

406 calculations in the stratosphere by RRTM (e.g., bands 6, 12-13 and 15-16). In 

407 contrast, MRAC incorporates the contribution of CO2 and H2O to the optical 

408 depth in these spectral bands. However, this contribution is usually rather 

409 small. 

410 Secondly, and more importantly, the differences between the stratospheric tem- 

411 perature profiles occur where RRTM has to use a temperature extrapolation 

412 for the absorption cross sections beyond the limits of its tabulated values. 

413 Fig. [5], shows the temperature profile for our test atmosphere Earth 1. Also 

414 shown is the validity range of RRTM as already indicated in Fig. [H This rep- 

415 resents the lower temperature limit for the tabulated absorption cross sections 

416 in RRTM, as stated above. 

417 As can be seen from Fig. El the calculated temperature values for the first test 

418 atmosphere are below the lower RRTM validity limit. Therefore, RRTM uses 

419 linear extrapolation to calculate the absorption cross sections. This introduces 

420 a large extrapolation error. On average, the calculated absorption cross sec- 

421 tions are a factor of 2-5 too low, depending on the spectral band. Sometimes, 
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422 the extrapolation performed by RRTM even yields negative absorption cross 

423 sections. The interpolation errors in MRAC, on the contrary, reach only 1-2% 

424 on average. 

425 Figure [6] shows the radiative fluxes and heating and cooling rates calculated by 

426 RRTM and MRAC in the test atmosphere Earth 1. Solar fluxes and heating 

427 rates differ by much less than 1 %. The thermal down- welling fluxes calcu- 

428 lated by RRTM and MRAC show large differences in the stratosphere below 

429 pressures of around 10~^-10~^ bar, reaching up to a factor of 5 in the upper 

430 stratosphere where pressures are below 10~^ bar. However, these differences 

431 are well below 1 W m~^, so are not discernible on the scale in Figure El The 

432 up-welling fluxes differ by only about 10 % in the stratosphere, since they 

433 are dominated by the tropospheric component. Hence, the calculated resulting 

434 cooling rates differ only by small amounts of 0.1-0.4 K day~^ and usually lie 

435 within 5-10 %, especially in the upper stratosphere. 

436 In order to assess the sensitivity of the model to errors in the absorption cross 

437 sections (hence, in optical depth and thermal fluxes), we artificially increased 

438 the optical depth in RRTM in the most important stratospheric band, the CO2 

439 15/im fundamental by factors of 2, 5, 10 and 20, respectively. Fig. [7] quantifies 

440 the effect of these sensitivity runs on the temperature profile. Fig. [7] a) shows 

441 the total optical depth calculated by RRTM and MRAC in the validation 

442 runs. Clearly, in the stratosphere, RRTM under-estimates the optical depth, 

443 as already discussed above. 

444 Fig. [7|b) shows the temperature profiles for the two Earth validation runs with 

445 MRAC and RRTM, as well as for a run with RRTM, but increased optical 

446 depth by a factor of 2. This factor of 2 is representative of the error in the 

447 absorption cross section calculations in the 15 //m band due to the required 

448 extrapolation in the T-p-range, as shown in Figure [H It can be seen that 

449 by increasing the optical depth in RRTM artificially, the temperature profile 

450 nearly reproduces the MRAC temperature profile. 
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451 These results show that conditions which differ too much from the Earth's 

452 standard atmosphere seem to pose problems for RRTM. This limitation was 

453 already noted in some of the previous studies performed with RRTM. Clearly, 

454 due to the use of an expanded temperature range, MRAC performs better than 

455 RRTM in these atmospheres. 



456 4 About the runs 



457 ^.i Absorption cross sections 



458 The absorption cross sections used in the runs performed for this work (sum- 

459 marized in Tables [3] and H]) were calculated assuming a N2-C02-background 



atmospher e, consisting of 95% mo 



cording to 



Kasting and Ackerman 



ecular nitrogen and 5% carbon dioxide. Ac- 



fll986h and 



TothI fl2000h . the foreign broad- 



462 ening coefficient for water is enhanced by a factor of 2 with respect to air for 

463 CO2 as a broadening gas and by a factor of 1.2 for N2 as a broadening gas. 

464 Similarly, for carbon dioxide the foreign broadening coefficient was enhanced 



465 by a factor of 1.3 when N2 was the broadening gas (IKasting and Ackermanl . 



19861 ). Accounting for the appropriate mixing ratios then yields an effective 



467 enhancement factor by which the foreign broadening parameter from HiTemp 

468 was multiplied before use in the cross section calculations described above. 

469 We compared the calculated cross sections of the assumed N2-C02-atmosphere 

470 (95% nitrogen and 5% carbon dioxide) with cross sections for some major spec- 

471 tral bands corresponding to different background atmospheres, i.e. different 

472 CO2 concentrations. These cross sections were obtained with the line-by-line 



radiative transfer model Mirart (ISchreier and Bottger . 



2OO3I ). The cross sec- 



474 tions from the 95%-N2-5%-C02-atmosphere agree within 5 % for most of the 

475 cases studied in this work, although for the runs with the lowest CO2 concen- 

476 trations, the agreement decreased to about 10 %. 
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477 



478 4-^ Model runs 



479 There were a total number of 12 nominal runs performed as shown in Table [31 

480 We assumed a constant background pressure of 0.77 bar N2 with variable 

481 amounts of CO2. Water vapour contents of the atmosphere were calculated as 
described in se ction 12.31 No other gases were present in the atmosphere, as in 



KastingI (119871 ). For several values of the solar constant {3=0.7, 0.75, 0.8, 0.85, 



484 0.9 and 0.95 present-day value), we increased the CO2 partial pressure until 



converged surface temperatures reached 273 K (runs 1-6) and 288 K 



runs 7- 



12), respectively. This procedure is similar to what was done by lKastingl (119871 ) 



The solar constant values were chosen as to loosely correspond to important 
events throughout the Earth's history, such as the beginning of the main se- 
quence life time of the Sun (5'=0.70, 4.6 Gy ago), the end of the late heavy 
bombardment (S'=0.75, 3.8 Gy ago), the rise of oxygenesis by cyanobacteria 
(iS'=0.80, 2.9 Gy ago), the first and the second oxidation event (S'=0.85 and 
>S'=0.90, 2 and 1.3 Gy ago respectively) and the Cambrian explosion (S'=0.95, 
0.6 Gy ago). Assigning geological ages to the solar constants used in the model 
(see Table [3]) is essential when comparing the calculated model CO2 concentra- 
tions to the available data and constraints on carbon dioxide. In this work, we 



496 chose t wo approximations of the so 



is from 



ar luminosity S{t) with time. The f irst one 



Caldeira and KastingI (119921 ) , the second approximation is from iGough 



498 ( 1l98ll ). The difference of these two formulations is most pronounced for earlier 



499 time periods before 1-2 Gy ago, although it is rarely larger than a few percent. 



500 The total surface pressure is always calculated from the relation psurf = PCO2 + 

501 + Pu20- Since we calculated surface temperatures of 273 K and 288 K, 

502 the amount of surface pressure which arose from the evaporation of water is 

503 about 6 mb and 17 mb, respectively. Note that the total surface pressure in 
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504 our model runs decreases with time, as we assume less CO2 partial pressure 

505 in our model atmospheres. There is no geological data available on the total 

506 surface pressure throughout time, however, our approach of constant partial 

507 press ure is consistent with assumptions made in previous studies (e.g. 



1987h . 



Kasting 



509 4-3 Parameter variations of surface albedo and relative humidity profile 

510 Previous works regarding the "faint young Sun problem" have concentrated 

511 on the increased greenhouse effect, as stated above in the Introduction. We 

512 also studied two important parameters affecting the surface temperature in our 

513 model, namely the surface albedo and the relative humidity (RH) distribution. 

514 The importance of cloud coverage for the surf ace temperature, e.g. on early 



515 Mars, has been studied by 



Mischna et al. 



(I2OO0I ). The effect in their model was 



516 quite large, yielding surface temperatures from 220 K to 290 K, depending 

517 on cloud optical depth and cloud location. In our model, the surface albedo 

518 simulates the presence of clouds, as stated above in the description of our 

519 model. 

520 Water vapour is a very effective greenhouse gas. In contrast to CO2, its content 

521 in the troposphere (where 99 % of the column resides) is controlled by the 

522 hydrological cycle (ocean reservoir, evaporation, subsequent condensation and 

523 precipitation) which is very sensitive to temperature. A critical parameter in 

524 climate models is the relative humidity parametrization, which clearly has a 

525 potentially large impact on the water content in the atmosphere and hence on 

526 the resulting greenhouse effect. 

527 In addition to the model runs described above, we therefore performed runs 

528 varying surface albedo and relative humidity profiles. These runs are based on 

529 runs 3 and 4 of Table 0, i.e. solar constants of S=0.8 and 15=0.85. Table H] 

530 summarizes the parameter variations. 
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531 In the runs from Table [3l the surface albedo is normally set to 74=0.21 (see 

532 Table Now, we assume surface albedos of A ±10 %, i.e. 0.19 and 0.23, 

533 respectively. 

culation of water vapour con - 



Manabe and Wetheraldl (119671 ) 



534 The relative humidity profile used for the ca 

535 centrations in Table [3] follows the approach of 

536 (referred to as RH=MW, see section 1231) . Here, for the additional parameter 

537 studies, we used three different RH profiles. The first one assumed a saturated 

538 troposphere (RH=1). The second profile used RH=0, i.e. water was removed 

539 from the atmosphere. These two profiles represent the lower and upper lim- 

540 its of the atmospheric water content. The third RH profile is a more realis- 



541 tic RH profile. It is based on a temperature correction to t 



he RH 



Manabe and Wetherald ( 



543 later used by 



19671) which wa s first proposed by 



profi le of 



Vardavas and Carverl (119851 ).: 



CessI (Il976h and 



R — R, 



surface 



_ ^l-0.03(Ts-288) 



(23) 



544 where /^surface • Rmw IS the relative humidity distribution from eq. 



This 



545 profile is referred to as RH=C. Figure M s. 



lows the difference between the 



546 two R H profiles from iManabe and Wetheraldl (I1967I ) and 



Vardavas and Carver 



547 ( 1l985l ). The total amount of wate r vapour is reduced by abou t 20 % by us- 

548 ing the temperature co r rectio n of IVardavas and Carverl (119851 ) compared to 



Manabe and Wetheraldl (Il967h . 



550 5 Results and Discussion 



551 5. 1 Examples of thermal structure and water profiles 



552 Figure [HI shows the temperature profiles for run 3 (plain line) and run 4 (dotted 

553 line). Run 3 (S* = 0.80) considered 20 mb CO2 partial pressure, run 4 (5 = 
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554 

555 
556 



0.85) roughly 3 mb (see Table [3]). These two runs were chosen because they 
represent interesting points in time, such as the advent of cyanobacteria (run 
3) and the first oxidation event (run 4). 



557 Clearly, in the troposphere the temperature differences for the two runs are 

558 rather small. In the lower to middle stratosphere (from around 10 to 25 km), 

559 run 4 (dashed line) shows lower temperatures than run 3 (plain line). In this 

560 region, absorption of solar radiation by CO2 is the dominant heating process, 

561 eventually causing a temperature inversion from the middle stratosphere up- 

562 wards (seen for both runs in Fig. [9] for higher altitudes above 25-30 km). In 

563 the upper stratosphere (above 25-30 km), where CO2 radiative cooling via the 

564 15 /im fundamental band sets in, run 3 shows lower temperatures than run 4. 

565 This reflects the lower flux and higher CO2 concentrations in run 3 compared 

566 to run 4. 



567 However, as already noted by ISegura et al.l (120031 ). the actual strength of the 



568 temperature inversion is not well determined due to less accurate transmittance 

569 calculations in the solar code and may well turn out to be a mere numerical 

570 problem rather than reflecting physical processes. 

571 Also indicated in Fig. [9] are the convective zones for both runs. The convective 

572 zone for the 5" = 0.8 case only reaches up to about 3 km altitude, whereas 

573 in the S = 0.85 case, it extends already to about 7.5 km, which is closer to 

574 the present-day, latitude-dependent value of 10-20 km. However, the lapse rate 

575 remains close to the adiabatic value for several kilometres after entering the 

576 regime of radiative equilibrium (above the dot-dashed lines in Fig. [H]). 

577 Figure [TU] is similar to Fig. [HI but it shows the resulting water proflles. There 

578 is considerably less water in the stratosphere for the S = 0.85 case due 

579 to enhanced condensation. The cold trap position is indicated, as well as 

580 the tropopause position, i.e. the boundary between the convective and non- 
581 convective regime. 
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5B2 Figures [9] and [lOl illustrate some interesting features regarding the thermal 

5B3 structure of the atmosphere. The cold trap, i.e. the highest point up to which 

5B4 water is allowed to condense out in the model, is no longer located near the 

585 tropopause, as is the case in the atmosphere of modern Earth. It is still as- 

5B6 sociated with a temperature inversion, but as this inversion occurs at higher 

5B7 altitudes, cold trap and tropopause (i.e., the boundary between convective and 

588 non-convective layers) are located at noticeably different altitudes. It seems 

589 that the observation that tropopause, cold trap and temperature inversion oc- 

590 cur at approximately the same height in the modern Earth's atmosphere is 

591 somewhat coincidental. 



592 5.2 Climatic constraints on CO2 partial pressures 



593 The results of the 12 nominal runs (see Table 4) are summarized in Fig.[TTl The 

594 lower plain line indicates the CO2 partial pressures corresponding to calculated 

595 surface temperatures of 273 K (runs 1-6), the upper plain line shows the partial 

596 pressures for surface temperatures of 288 K (runs 7-1 2). For compar ison, the 

597 dashed line shows the partial pressures as calculated by lKastingl (119871 ). Fig. [TT] 

598 shows that our new radiative scheme requires considerably le ss CO^ (abou t a 



599 factor of 2-15) to achieve an ice- free surface than the study of iKastind (119871 ) 



600 Note that the radiative transfer scheme used by lKasting (119871 ) was applicable 

601 to this type of atmospheres, as is the one we used. 



602 Also indicated in Fig. [TT] are the upper limits of CO2 partial pressures as 

603 inferred from sedimentary data. 



604 Several attempts have been made to constrai n the CO9 conte n t of the early 



Earth atm osphere. For 



example, among others. 



Rye et al 



(119951) 



Hessler et al 



606 ( I2OO4J ) and lTowd (1l985l ) have placed upper limits on CO2 partial pressures dur- 

607 ing different periods of the mid- Archaean to early Proterozoic of approximately 

608 10-100 PAL (Present Atmospheric Level), depending on temperature. These 
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609 limits we re based on exp erime ntal sediment data and on biological consid- 



610 erations. 



Rye et al. 



(119951) and 



Hessler et al. 



(120041 ) argued for upper limits 



on CO2 partial pressures ba sed ori the o 



Dservations that in ancient Archaean 



612 rocks, siderites are missing. 



Towd (119851 ) stated atmospheric upper levels of 



613 CO2 based on the observation that anaerobic photosynthesis and nitrogen fix- 

614 ation (which are both believed to be evolutionary ancient) are incompatible 

615 with high CO2 concentrations. 



616 Previous climate studies of the early Earth's atmosphere calculated very high 

617 CO2 partial pressures in order to raise the surface temperature above 273 K. 

618 The model CO2 values were generally above 50 mb well into Proterozoic age 

619 (e.g., Kasting 1987, see also Fig. [TTl) . However, the experimental data from 

620 sediments and biology (see above) were of the order of several mb, i.e. an 

621 order of magnitude lower. This contradiction is known as the " faint young Sun 

622 problem" . 



623 Fig. [TT] however shows that for the late Archaean (solar constant S=0.85), our 

624 new results are compatible with the paleosol records, hence the "faint young 

625 Sun problem" might be resolved for this time period. 



626 One should note that the absence of methane, ozone, ammonia or other green- 

627 house gases in our results implies that we calculated only a lower limit for 

628 surface temperatures. Photochemical models of the anoxic Archaean and low- 

629 oxygen Proterozoic atmospheres ha ve shown tha t meth ane could build up to 



630 concentrations of the order of 10 ^ ( 

631 contributed to the greenhouse effect. 



'avlov et al 



Haqq-Misra et al. 



20031 ) whic h could also have 



(I2OO8I ) proposed higher 



632 hydrocarbons such as C2H6 in a methane-rich atmosphere as radiative gases. 

633 Those effects were not investigated here. Our studies imply, however, that 

634 much less to none additional greenhouse gases are required to warm the early 

635 Earth. 
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636 5.2.1 Effect of parameter variations of surface albedo and relative humidity 



637 Table m summarizes the results of the parameter variations. Shown are the val- 

638 ues of CO2 partial pressure required to obtain the desired surface temperature 

639 of 273 K. 

640 Upon lowering the surface albedo value from 0.23 to 0.19, the partial pressure 

641 of CO2 required to keep the surface at 273 K is reduced from 28.5 mb to 11.0 

642 mb (at 5=0.80) and from 5.5 mb to 1.5 mb (at >S'=0.85). This is a lowering of 

643 the amount of CO2 by factors of 2.6 and 3.7, respectively. 

644 When changing the RH profile from saturated (RH=1) to water-free condi- 

645 tions (RH=0), the necessary partial pressure of CO2 has to be increased from 

646 4.9 mb to 266.8 mb (5=0.80) and from 0.6 mb to 180.6 mb (5=0.85). How- 

647 ev er, like expected, the compar ison betwe e n the more realistic RH profiles 



of 



Manabe and Wetheraldl (119671 ) and ICessI (119761 ) yields much smaller differ 



649 ences. When changing the RH profile from RH=MW to RH=C, CO2 partial 

650 pressures must be increased from 19.1 mb to 27.3 mb (5=0.80) and from 2.9 

651 mb to 4.9 mb (5=0.85). These differences amount to 43 % and 69 %, respec- 

652 tively, which is much lower than the differences obtained by varying the surface 

653 albedo. 



654 These parameter variations demonstrated that the surface albedo, hence clouds, 

655 is a very important parameter in view of the surface temperature. The results 

656 obtained imply that the latter is probably more important than the RH profile, 

657 although the RH profile has to be incorporated more consistently to accurately 

658 constrain CO2 partial pressures for the "faint young Sun" problem. This calls 

659 for a more elabor ated ID model incorporatin g clouds and cloud formation as 



660 was done by, e.g. 



Colaprete and Toon! (120031 ) for early Mars. In the future, we 



661 plan to add clouds to our code. 

662 The results of the parameter variations, however, did not significantly change 

663 our main conclusion, namely that the CO2 values for the late Archaean cal- 
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664 culated by our improved model are consistent with observational data. The 

665 CO2 values obtained for the parameter variations range between 1.5-5.5 mb 

666 of partial pressure. This is still close to the values inferred from the paleosol 

667 record. 

668 6 Summary 

669 In this work, we addressed the "faint young Sun problem" of the ice-free early 

670 Earth. In order to do this, we applied a one-dimensional radiative-convective 

671 model to the atmosphere of the early Earth. Our model included updated 

672 absorption coefficients in the thermal radiative transfer scheme. 

673 The validations done for the new radiative transfer scheme have been described. 

674 The new scheme is found to perform significantly better than a previous scheme 

675 under conditions deviating from the modern Earth atmosphere. 

676 We then studied the effect of enhanced carbon dioxide concentrations and 

677 parameter variations of the surface albedo and the relative humidity profiles 

678 on the surface temperature of the early Earth with the improved model. 

679 Our new model simulations suggest that the amount of CO2 needed to keep 

680 the surface of the early Earth from freezing is significantly less (up to an order 

681 of magnitude) than previously thought (see Figure fTTil . 

682 For the late Archaean and early Proterozoic period (around 2-2.5 Gy ago), the 

683 calculated amount of CO2 (2.9 mb partial pressure) which is needed to obtain 

684 a surface temperature of 273 K is compatible with the amount inferred from 

685 geological data, contrary to previous studies (see Figure [TTl) . The apparent 

686 contradiction between model constraints on CO2 and sediment data disappears 

687 for this time period. 

688 Upon varying model parameters such as the surface albedo and relative hu- 
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6B9 midity profile, we found tliis conclusion to be robust. The calculated CO2 

690 partial pressures for the late Archaean (1.5-5.5 mb) are still consistent with 

691 the geological evidence. 
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837 Figure captions 



838 Figure [T] 

839 Range of T-p values used to obtain absorption cross sections in the two radia- 

840 tive schemes RRTM (hght grey) and MRAC (dark grey). 



842 Figure E) 



843 Cumulative k distribution for parts of the 6.3 urn. H2O funda r nenta l band 



844 T=240 K and p=10 mb: Comparison between 

845 line) and our algorithm (dotted). 



Lacis and OinasI (Il99ll ) (plain 



847 Figure [3l 

848 Cumulative k distribution for parts of the 15 /xm CO2 fundamental band 

849 between 630 a nd 70 cm~^, T=260 K and p=507 mb: Comparison between 



Mlawer et al. 



(119971 ) (plain line) and our algorithm (dotted). 



852 Figure HI 

853 Comparison of validation temperature profiles (RRTM plain line, MRAC dot- 

854 ted hne). 



856 Figure El 

857 Limits of the RRTM temperature grid. Atmospheric conditions as for the stan- 

858 dard CO2 case in fig. H RRTM profile (plain), MRAC (dotted). The shaded 

859 area designates the validity range of RRTM as already indicated in fig. [TJ 
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861 Figure E) 



862 Flux profiles (thermal and solar up- and down-welling fluxes) (left panel) and 

863 heating and cooling rates (right panel) calculated by the two radiative schemes 

864 in the test atmosphere. 

865 

866 Figure [T) 

867 Effect of variations in the optical depth on the temperature proflle calculated 

868 by RRTM: Left panel, total optical depth, right panel, temperature proflles. 

869 

870 Figure [HI 

871 Differences of the two relative humidity proflles used for the parameter varia- 

872 tions 

873 

874 Figure M 

875 Temperature proflles for runs with solar constants of = 0.8 (run 3, plain 

876 line) and S = 0.85 (run 4, dotted); the background pressure is 0.77 bar of N2. 

877 The convective zones are indicated for both runs. 

878 

879 Figure [101 

880 Water proflles for runs with solar constants of S* = 0.8 (run 3, plain line) and 

881 S* = 0.85 (run 4, dashed); the background pressure is 0.77 bar of N2. Cold trap 

882 positions and tropopause positions are indicated for both runs. 

883 

884 Figure [TT] 

885 Minimal values of CO2 partial pressure required to obtain chosen surface tem- 
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8B6 peratures of 273 K and 288 K at a fixed N2 partial pressure for different solar 



constants: shown are the curves for the new model (p 



288 K, lower line 273 K) and the model of lKastingI (119871 ) (dashed)). Symbols 



ain lines, upper line 



889 are included which represent the upper CO2 limits derive d from the sediment 

890 rec ord age c onversion by ICaldeira and Kastingj (119921 ). A: age conversion 

891 bv boughl Jl98lf ). 
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892 Tables 



Table 1 

Contribution of model species to the temperature profile via radiative transfer for 
solar or thermal radiation, adiabatic lapse rate formulations or heat capacity con- 
tributions (x: active species, -: inactive species) 
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Table 2 

Initial values, boundary conditions and parameters in the model (IV: initial value, 
BC: boundary condition, P: parameter) 



Quantity 


Value 


Type 


To-profile 


US Standard 1976 


IV 


TOA incident flux 


solar spectrum 


BC 


TOA po 


6.6-10-5 bar 
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Surface albedo 
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Zenith angle 
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P 
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Table 3 

Summary of model runs performed for this work (partial pressures p^^ PCO2 
and surface pressure Psurf ^-re in mb, surface temperature Ts^rf in K) 
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Table 4 



Summ ary of sensitivi ty runs performed (RH profiles: MW iManabe and Wetherald 
11963), C [c^(fl976l )l. -•artial CO2 pressure PCO2 reach 273 K in mb. 



Run number 
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